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Quantitative polymerase chain reaction-high-resolution melting (qPCR-HRM) analysis was used to screen for 
mutations related to drug resistance in Mycobacterium tuberculosis. We detected the C526T and C531T mutations 
in the rifampicin resistance-determining region (RRDR) of the rpoB gene with qPCR-HRM using plasmid-based 
controls. A segment of the RRDR region from M. tuberculosis H37Rv and from strains carrying C531T or C526T 
mutations in the rpoB were cloned into pGEM-T vector and these vectors were used as controls in the qPCR-HRM 
analysis of 54 M. tuberculosis strains. The results were confirmed by DNA sequencing and showed that recombinant 
plasmids can replace genomic DNA as controls in the qPCR-HRM assay. Plasmids can be handled outside of bio- 
safety level 3 facilities, reducing the risk of contamination and the cost of the assay. Plasmids have a high stability, 
are normally maintained in Escherichia coli and can be extracted in large amounts. 
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The development of rapid and affordable assays for 
diagnosing multidrug-resistant tuberculosis (TB) is 
necessary to control the spread of the disease (Bwanga 
et al. 2009). The rapid detection of Mycobacterium tu- 
berculosis reduces treatment failure, the rate of trans- 
mission and the period of infectiousness (Martin et al. 
2008, Koh et al. 2012). 

High-resolution melting (HRM) analysis can be used 
to distinguish DNA sequences with slight differences 
in the melting curve profile. This technique measures 
the fluorescence of a DNA dye released from double- 
stranded DNA with high accuracy (Ghorashi et al. 2010). 
HRM analysis has been used in screens for gene muta- 
tions prior to DNA sequencing, in methylation studies 
and for genotyping, with high throughput, reliability, ra- 
pidity, simplicity and cost efficiency (Ereqat et al. 2010, 
Kagkli et al. 2012). 

A number of DNA deletions, insertions and single- 
nucleotide polymorphisms in the M. tubercidosis ge- 
nome are well established as causing resistance to anti- 
TB drugs (Nimri et al. 2011, Malik et al. 2012). HRM 
analysis of these mutations are useful for the rapid pre- 
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diction of the drug sensitivity profile and, consequently, 
allow for the adjustment of treatment regimens (Ramirez 
et al. 2010, Lee et al. 2012, Yadav et al. 2012). 

Rifampicin is considered a good marker for multi- 
drug resistance. This drug binds to the P-subunit of RNA 
polymerase encoded by rpoB and inhibits DNA-directed 
RNA synthesis. Mutations in the rpoB gene are associ- 
ated with rifampicin resistance and are usually located 
within an 81 bp hypervariable region of the rpoB gene 
known as the rifampicin resistance-determining region 
(RRDR) (Khadka et al. 2011, Nimri et al. 2011). Quan- 
titative polymerase chain reaction-HRM (qPCR-HRM) 
analysis of the RRDR locus can be used to rapidly screen 
for rpoB mutations (Ong et al. 2010, Yadav et al. 2012). 

To improve the specificity of HRM analysis, ge- 
nomic DNA samples with specific mutations are used 
as positive or negative controls (Choi et al. 2010, Ong 
et al. 2010). However, the production of large amounts 
of genomic DNA carrying resistance-associated muta- 
tions must be performed in a biosafety level 3 laboratory. 
Cloning strategies that allow for the large-scale produc- 
tion of recombinant DNA in Escherichia coli or other 
less pathogenic bacteria may provide a useful method to 
obtain M. tubercidosis mutant sequences that will im- 
prove the performance of HRM analysis. 

In this study, plasmids carrying the C526T or C531T 
mutation or no mutations in the RRDR region of rpoB 
were cloned into the pGEM-T easy vector and used 
as positive and negative controls, respectively, in the 
HRM assays. 

A total of 54 M. tuberculosis isolates were obtained 
from our M. tuberculosis collection. Rifampicin sus- 
ceptibility was determined by the resazurin microtitre 
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assay (REMA) plate method according to Palomino et 
al. (2002). The rifampicin susceptibility profile was also 
obtained by the nonradiometric BACTEC MGIT 960 
method (BD, USA) (Miyata et al. 2011). 

Mycobacterial DNA was extracted using a QIAamp 
DNA mini kit (Qiagen, USA) according to the manufac- 
turer's instructions. The primers RpFw 5'-GTCGGC- 
GAGCTGATCCAAAAC-3' and RpRv 5'-GGTACG- 
GCGTTTCGATGAACC-3' (IDT, USA) were designed 
using Primer Premier software 5.0 (Biosoft, USA) to 
amplify a 409-bp segment of rpoB containing the RRDR 
region. PGR assays were performed with a PTC 200 
thermocycler (MJ Research, USA) under the following 
cycling conditions: (i) an initial step at 95°C for 10 min, 
(ii) 30 cycles at 95°C for 45 s, 58°C for 30 s and 72°C for 
1 min and (iii) a 10-min step at 72°C. 

The PGR products were purified with a QIAquick 
PGR purification kit (Qiagen, USA) and sequenced with 
an ABI 3730 DNA Analyser using the BigDye Termi- 
nator v3.1 Gycle Sequencing Kit (Applied Biosystems, 
USA) and the previously described set of primers. The 
sequence was analysed with BioEdit software v7.1.3 (Ibis 
Biosciences, USA). The 409-bp segment was cloned into 
the pGEM-T Easy Vector (Promega, USA) and plasmid 
DNA was extracted with the Plasmid Midi kit (Qiagen, 
USA) according to the manufacturer's protocol. 

To perform qPGR-HRM, the primers HRM-Fw 5'- 
TGGGGGGGATGAAGGAGTT-3' and HRM-Rv 5'-GT- 
GGAGGTGGGGGAGGTG-3' (IDT, USA) were designed 
to amplify a 158 bp segment oirpoB spanning the RRDR 
region. qPGR-HRM was performed in duplicate on a Ro- 
tor-gene 6000 platform (Qiagen, USA) in a reaction mix- 
ture containing 10 |J,L of 2X Amplitaq Gold PGR Master 
Mix (Applied Biosystems, USA), 1 |.iL of DNA, 200 nM 
each primer, 0.6 |iM Syto9 (Invitrogen, USA) and PGR- 
grade water to a reaction volume of 20 (xL. 

Recombinant plasmids and M. tuberculosis DNA 
were diluted to obtain qPGR threshold cycle (Gt) values 
between 19-29 (Ramirez et al. 2010). The cycling param- 
eters were as follows: one step of 10 min at 95°G and 40 
cycles of 20 s at 95''G and 1 min at 62°G. HRM was per- 
formed with an initial step at 60°G for 1 min followed by 
temperature increases over the range of 65-95°G at a rate 
of 0.1°G/s with constant fluorescence acquisition. 

The PGR and HRM results were analysed using Ro- 
tor-gene 6000 software v. 2.0.2 (Qiagen, USA). The melt- 
ing curves were analysed using Rotor-gene ScreenGlust 
HRM software (Qiagen, USA). This software was used 
to optimise the assay and group the samples. The DNA 
sequencing of the recombinant plasmids pGEM-G526T, 
pGEM-G531T and pGEM-wild type (WT) containing the 
rpoB coding segment confirmed the expected sequences 
and no additional mutations that could interfere with the 
melting curve profiles were found. 

The concentrations of plasmid and genomic DNA 
with Gts between 19-29 provided accurate data for the 
qPGR-HRM analysis. Similar results were obtained by 
Ramirez et al. (2010) using genomic DNA. PGR products 
that were directly amplified from M. tuberculosis H37Rv 
genomic DNA presented the same melting curve profile 
as pGEM-WT. Moreover, it was possible to discriminate 
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Fig. 1: quantitative polymerase chain reaction-high-resolution melt- 
ing profiles of recombinant plasmids. Normalized fluorescence vs. 
temperature of the 158 bp segment from Mycobacterium tuberculosis 
genomic DNA and recombinant plasmids. 1: pGEM-C526T and ge- 
nomic DNA of M. tuberculosis carrying mutation C526T; 2: pGEM- 
WT and M. tuberculosis U^^K genomic DNA; 3: pGEM-C531T and 
genomic DNA carrying mutation C531T. 
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Fig. 2: quantitative polymerase chain reaction-high-resolution melting 
profiles of DNA samples. Differential plots of the 158 bp segment of 
rpoB from Mycobacterium tuberculosis genomic DNA and recombi- 
nant plasmids. 1 : pGEM-wild type and M. tuberculosis H^^R^ genomic 
DNA, C526T; 2: pGEM-C526T and genomic DNA of M. tuberculosis 
samples carrying mutation C526T, C531T; 3: pGEM-C531T and ge- 
nomic DNA samples carrying mutation C531T. Variations: samples 
presenting other mutations. 



pGEM-WT from pGEM-G526T, pGEM-G531T and M. 
tuberculosis genomic DNA with other mutations. In 
addition, the recombinant plasmids pGEM-G526T and 
pGEM-G531T had the same melting curve profiles as 
M. tuberculosis genomic DNA carrying the G526T and 
G531T mutation, respectively (Fig. 1). 

Visual discrimination of the melting curve profiles of 
pGEM-G526T and pGEM-G531T was difficuh due to the 
nucleotide change being the same and the proximity of 
the mutations in the DNA strand. However, the software 
used to analyse data was able to accurately distinguish 
these mutations. The detection of slight differences in 
the melting curve profile of two different mutations is 
important to specifically detect only mutations associ- 
ated with phenotypic drug resistance profiles (Figs 1, 2). 
These results revealed that recombinant plasmids can 
replace genomic DNA as positive and negative controls 
in the detection of the two main mutations in the RRDR 
region correlated with phenotypic rifampicin resistance 
without interfering with the data interpretation. 
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TABLE 



Rifampicin phenotypic profile of Mycobacterium tuberculosis strains and detection 
of mutations (C526T and C531T) in the rifampicin resistance-determining region obtained 
by quantitative polymerase chain reaction-high-resolution melting (qPCR-HRM) and confirmed by DNA sequencing 



Sequence 


qPCR-HRM 


DNA 
sequencing 


Rifampicin phenotypic profile 
(resazurin microtitre assay ) 

Susceptible Resistant 


Mutation C526T 


3 


3 


3 


Mutation C53 IT 


15 


15 


15 


Other mutations 


12 


12 


2 10 


No mutation 


24 


24 


18 6 


Total of strains 


54 


54 


54 



The plasmid-based controls were used to screen for 
C526T and C531T mutations in the RRDR region of 54 
M. tuberculosis strains. The qPCR-HRM results gave 
four curve profiles that allowed for the mycobacterial 
strains to be divided into four groups: (i) WT, (ii) C531T 
mutation, (iii) C526T mutation and (iv) other mutations. 
This fourth group contained all of the samples with mu- 
tations other than C526T and C531T (Table). 

The RRDR regions of all strains were sequenced. 
The results confirmed that 100% of the strains with mu- 
tations were correctly typed as carrying the C526T or 
C531T mutation or other mutations (Table). The plasmid- 
based controls had high specificities; pGEMT-WT was 
able to differentiate all WT strains from strains with nu- 
cleotide changes. Strains presenting the C526T or C531T 
mutation were easily identified using their respective 
plasmids and no strains with other mutations were in- 
cluded in groups 1, 2 or 3. 

Mutations in the RRDR region are usually found in 
clinical isolates resistant to rifampicin and most of these 
mutations occur in codons 516, 526 and 531 (Lingala et 
al. 2010, Khadka et al. 2011, Yadav et al. 2012). In this 
study, six strains with no mutations in the RRDR region 
of rpoB, but that were rifampicin resistant according to 
the REMA, were identified as WT by the HRM analysis 
and DNA sequencing. This resistance could be due to 
mutations in other regions of the M. tuberculosis genome 
that were not evaluated in this study (Siu et al. 2011). In 
addition, two strains presenting mutations in the RRDR 
and included in group 4 were susceptible to rifampicin 
according to the REMA, but resistant to rifampicin ac- 
cording to the nonradiometric BACTEC MGIT 960 gold 
standard method (Miyata et al. 2011). DNA sequencing 
showed a single mutation (C531G) in the RRDR of one 
strain, whereas two mutations (C505T and C526T) were 
found in the second strain. These findings confirm the 
BACTEC MGIT susceptibility results and show that the 
C531G and C526T mutations are related to rifampicin 
resistance (Mark et al. 2005). Additionally, the presence 
of the C531G or C505T mutation in the 158 bp segment 
generated a specific melting curve profile for each strain 
and, therefore, these strains were correctly classified 



into group 4 by the qPCR-HRM analysis. The inclusion 
of new plasmid controls for other nucleotide changes 
that cause rifampicin resistance will certainly reduce the 
number of samples placed in group 4. 

In contrast to some molecular methods, HRM analy- 
sis can be easily adapted to accurately detect mutations 
according to the necessities of each region or country. 
Plasmid libraries can be developed with gene mutations 
that are frequently related to drug resistance and stored 
for a long time in large quantities. These libraries can 
be easily maintained in E. coli, which can be grown in 
inexpensive media. There is no need to frequently grow 
drug-resistant M. tuberculosis strains exclusively for ex- 
tracting DNA samples for use as reference sequences. 

The detection of drug resistance with qPCR-HRM 
using plasmid controls requires less than 3 h for comple- 
tion, whereas phenotypic assays usually take more than 
a week. A rapid test to detect drug resistance has the po- 
tential to greatly improve the quality of TB treatment be- 
cause the patient will receive the correct drug, the treat- 
ment cost will be lower and the probability of spreading 
multidrug-resistant TB will be reduced. 

In conclusion, plasmid-based controls are a good al- 
ternative for the routine detection of known mutations 
associated with drug resistance using qPCR-HRM. 
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